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Using the theory of diffusion in graphs, we propose a model to study mesoscopic transport through a dif- 
fusive quantum dot. The graph consists of three quasi-lD regions: a central region describing the dot, and 
two identical left- and right- wires connected to leads, which mimic contacts of a real system. We find the 
exact solution of the diffusion equation for this graph and evaluate the conductance including quantum cor- 
rections. Our model is complementary to the RMT-models describing quantum dots. Firstly, it reproduces 
the universal limit at zero temperature. But the main advantage compared to RMT-models is that it allows 
one to take into account interaction-induced dephasing at finite temperatures. Besides, the crossovers from 
open to almost closed quantum dots and between different regimes of dephasing can be described within a 
single framework. We present results for the temperature dependence of the weak localization correction to 
the conductance for the experimentally relevant parameter range and discuss the possibility to observe the 
elusive OD-regime of dephasing in different mesoscopic systems. 
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1 Introduction 

In the last decades, dephasing in quantum dots has been studied experimentally and theoretically in great 
detail. The theoretical description is largely based on results from random matrix theory (RMT), empha- 
sizing the universality in the description of a dot, when spatial degrees of freedom become negligible. 
While the universal limits are well understood and reproduced in many experiments, a prediction of the 
full temperature dependence of quantities which are sensitive to dephasing, such as quantum corrections 
to the classical conductance, Ag, are challenging existing theories. Since RMT is not able to describe the 
T dependence on its own, several extensions were introduced in the past to describe their dependence on a 
dephasing time r^, which has to be included phenomenologically, see Section 2 for details. 

One of the well-know problems in the theory of dephasing in quantum dots originated from the pre- 
dictions of a seminal paper by Sivan, Imry and Aronov, who showed that dephasing in the so-called OD 
regime {T <C i?Th, where iSxh is the Thouless energy), behaves as ~ T^^, which results from Pauli 
blocking of the Fermi sea [1]. However fundamental the origin of OD dephasing is, it has so far not been 
observed experimentally. One possible reason for this might be the fact that dephasing is very weak in this 
regime, such that quantum corrections may reach their universal Umit Ay ^ 1. In general, if the dephasing 
time is much larger than the time the electron spends in the dot, Ag is governed by a dwelling time Tdw 
and becomes almost T independent. The remaining small T-dependent part of Ag can be masked, for 
example, by other T-dependent effects coming from contacts or leads. Thus, to facilitate an experimental 
observation of OD dephasing, a comprehensive theory of transport in the quantum dot connected to leads 
via some contacts is needed, which goes beyond the simple picture provided by RMT. 

In this paper we propose an alternative to the RMT description of the quantum dots. Namely, we foUow 
the ideas of [2, 3] and model the quantum dot as a network of ID wires and use the theory of diffusion 
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in graphs to calculate r^p and Ag. Earlier papers either focused only on small graphs, such as ID rings 
[4, 5, 6], or the authors introduced only phenomenologically [2, 3, 7]. We generalize the theory of 
for arbitrary graphs and include the regime T < i?Th by taking into account the Pauli principle. Using this 
theory, we calculate for a network describing a quantum dot, taking into account effects of the contacts 
and the leads. This allows us to demonstrate that the T^-dependence of the dephasing rate in OD regime is 
substantially distorted in usual transport measurements in quantum dots. 

The rest of the paper is organized as follows: In Section 2 we give a brief review of known results for 
dephasing in quantum dots. In Section 3, basic results from the theory of diffusion in graphs are presented, 
and in Section 4 we will apply this theory to construct a solvable quantum dot model as an alternative to 
the well-known RMT models. Results for the quantum corrections to the conductance and the dephasing 
time are presented in the following Sections. In the conclusions we compare different experimental setups 
where OD dephasing could be observed. 



2 Dephasing in quantum Dots: Brief review of known results 

It is well-known that the conductance g of a disordered normal metal is reduced due to quantum mechanical 
interference of the electron wave functions scattered at static impurities. It has been found that the reduc- 
tion of g can be expressed via the return probability of coherent electron paths, P{x, x, t), (the so-called 
Cooperon) integrated over time and space [8]: 

Ag = g - go = -4ETh dt J d'^x P{x,x,t) . (1) 

Here go is the classical conductance measured in units of e'^/h , Etu = D/Vl^ is the Thouless energy of 
the system, D is the diffusion constant and Q, is the largest size of the system. Ag is usually referred to as 
the weak localization correction. 

Quantum coherence is suppressed by a constant magnetic field and by time-dependent (noisy) fields, 
or when closed electron paths contributing to P{x, x, t) in Eq.(l) are dephased due to inelastic scatter- 
ing events. The time-scale associated with the latter is called dephasing time t^. In the absence of other 
sources of dephasing, yields an infrared cutoff for the time-integral, Eq.(l), and governs the temper- 
ature dependence of A^ [9]. At low temperatures, T < IK, where phonons are frozen, is dominated 
by electron interactions and depends on the dimensionality d and the geometry of the system. The T- 
dependence of in different regimes is governed by an interplay of with the thermal time tt = 1/T 
and the Thouless time rxh = 1 /£^Th, see Table 1 for a summary of known regimes in ID and 2D [10]. For 



Table 1 Dephasing rate as a function of temperature T. 
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low temperatures and small system sizes, when Et\i is the largest energy scale, dephasing becomes effec- 
tively zero-dimensional (OD). Therefore, it must be relevant for transport in metallic (diffusive or chaotic) 
quantum dots [ 1 ] . 

Note that OD dephasing requires confinement of the electron paths during times larger than ttii, since 
quantum corrections become T independent for T^p ^ rxh in fully open systems. As an example, consider 
the case of a quasi- ID wire of length L connected to absorbing leads, where Ag reads [11]: 



A,g = -4 



(2) 
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Fig. 1 (a) A graph consisting of 9 wires and 6 vertices, denoted by Greek letters, (b) A quantum dot realized as a 
graph with dimensionless coordinate Z = xjL. The labels Zi denote the position of the leads (i = 0,9) and the 
vertices (i = 1 ... 8) on the scale Z. Furthermore, the numbers i = 1 . . . 8 correspond to the ith row or column of the 
vertex matrix Eq.(lO) 



Thus, a detailed calculation of Ag including requires solving the full diffusion equation of the connected 
quantum dot, which is hard to achieve analytically for confined systems. 

One way to circumvent this problem is to apply random-matrix theory (RMT) to the scattering matrix 5, 
describing transmission and reflection in the sample. In such an RMT-model one assumes that the elements 
of the Hamiltonian H describing the systems are either real (Gaussian orthogonal ensemble, /?goe = 1) 
or complex (Gaussian unitary ensemble, /3gue = 2) random numbers corresponding to a system with 
time-reversal symmetry or broken time-reversal symmetry. ' Imposing a Gaussian probability distribution 
P{H), the scattering matrix S can be constructed using so-called R-Matrix theory. Alternatively, a sim- 
pler approach starts from a probability distribution of the scattering matrix directly, which is of the form 
P{S) = const, and S is again only restricted by symmetry arguments. From the scattering matrix, the 
full non-pertubative distribution of the transmission matrix and the conductance can be obtained. While 
RMT is unable to predict the temperature dependence of Ag on its own, the difference in g of the cases 
Pgoe and /3gue is equivalent to Ag in the limit of T 0. Extensions to RMT have been introduced in 
the past to describe the dependence of Ag on a dephasing time [12], e.g. by including a fictitious voltage 
probe into the scattering matrix which removes electrons from the phase-coherent motion of the electrons 
in the quantum dot [13], or by including an imaginary potential equal to —i/2T^ in the Hamiltonian from 
which the scattering matrix is derived [14]. It is expected that included in such an approach has the 
same form as stated in Table 1 for T <C i?Th, i-C- V oc T^^, but a proof of this expectation and a theory 
of a crossover between different regimes is still missing. 



3 Diffusion in graphs 

In this section, we present basic results from the theory of diffusion in graphs, following Ref.[l 1]. A graph 
is defined as a set of quasi-lD wires connected to each other at vertices, see the example shown in Fig. 1(a). 
In this section we will show how the solution to the Laplace transformed diffusion equation, 

{^-DA)P^{x,v)^5{x~y), (3) 

between arbitrary vertices (with coordinates x and y) of such a graph can be obtained. The time-dependent 
probability, required to calculate Ag and Tip, can be obtained via an inverse Laplace transform: 

-j p+ica 

P{x,y,t) = — d7e^*P^(x,y). (4) 

It is convenient to introduce the following quantities: We denote the wire between arbitrary vertices a 
and /3 as {a/3) and its length as Lap. Furthermore, the running coordinate along this wire (measured from 
a) is denoted Xa/s, and in the following, we will not distinguish a vertex from the coordinate of the vertex 
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on the graph: For example, P{a, y) is equivalent to Mmx^^^^o P{xai3,y), for any neighboring vertex /3 of 
a. The current conservation at some vertex a can be written as follows: 

where the symbol means summation over all wires (a/3) which are connected to a. 

Consider the point x lying at the coordinate Xap of wire (a/3) in Fig. 1(a). The probability to reach 
X from some arbitrary other point y of the graph can be expressed in terms of the probabilities from the 
neighboring vertices of x, i.e. a and /3: 

Pj {y, a)sinh ( y^j/D {Lap - Xap)) + Pj{y, P)smh( ^/J/D Xaf}) 

Pf (y, ^) = ^ , , _ X ^ ^ • (6) 



smhiy/j/D La 



Validity of the solution (6) can be checked directly by substituting Eq.(6) into Eq.(3). 
Inserting (6) into (5) yields the following equations for vertex a: 

P^iti,a)Y,Vl7Dcoth(y^Laf3) - Y.P^{^i,P) 7^^^ ^=D5a^f.- (7) 

(a/3) (a/3) slnh ( ^7/DL„/3 ] 



Writing down Eq.(7), for every vertex of the graph, we obtain a set of linear equations which can be solved 
for arbitrary vertices. Let us define a matrix M'' as follows: 

Mlp = (SapVl/Dcoth(^./^ Las) - S5fiVl/Dsmh(^^/^Las) j . (8) 

It is easy to check that the diffusion probability between arbitrary vertices of the graph is given by the 
entries of the inverse matrix divided by the diffusion constant [7, 11]: 



^7(«>/3) = ;^(-^^U- (9) 



4 A graph model for a connected quantum dot 

In this section we explain how to describe a connected quantum dot by a network of ID wires. The main 
advantage of this model is that an exact solution to the diffusion equation can be found. 

Consider the network shown in Fig. 1(b). It includes 8 vertices and describes a quantum dot of total 
length 2L attached via two contacts of length L to absorbing leads. Multiple wires connecting the same 
vertices (e.g. the three wires connecting vertex 4 with vertex 5) mimic a larger number of channels. Below, 
we use a dimensionless coordinate Z = xjL; the position of the leads is fixed at Zo = —2, Zo, = +2 
and the position of the 3rd and 6th vertex, describing the connection of the dot to the contacts, is fixed at 
Z-j, = —1, Zg = +1. The remaining 6 vertices are auxiliary: There are 3 regions in the system marked by 
"L" (left contact), "D" (dot) and "R" (right contact). We would like to describe diffusion from an arbitrary 
point in the system to another. Therefore, we have to place two additional vertices in each region L, D, R. 
Positions of these vertices define running coordinates. They are arbitrary within the corresponding region, 
thus each region is subdivided into 3 wires of varying length. The running coordinates can be expressed 
via the length of the connecting wires, e.g. the length of the wire connecting vertices 1 and 2 is given by 

(Z2 - Zi). 

To describe confinement of the electrons, we assume that all vertices in the regions L and R (including 
boundaries) are connected by single wires while the vertices in the dot (including its boundaries) are 
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connected via Na wires. This allows us to tune the system from a simple wire at Nd = 1 to an almost 
closed quantum dot for Nd oo. The corresponding vertex matrix M'' , defined in Eqs.(8), is given by 
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(10) 



We have introduced 2x2 blocks. 



Ml - Ai;ji23 , 



Ml 



and 



Ml = ^1^789 . 



(11) 



which are given by: 



Ml,i 



C0th(V7[^J- - Z^]) +COth(\^[Zfe - Zj] 

-l/sinh(V7[^fc-^j]) 



-l/sinh(v^[Zfc-Z,]) 



(12) 



Expressions for the entries "S"' and "C", which correspond to connected vertices, read 

Sl,r = -l/sinh(v^(-l T ^2;7)) , Sld;DR = -iVd/sinh(V^(l ± Z^.^^ , 

ClD;DR = NdCOth(^{l ± Z4;5)) + COthf v/^(-l T ^2;7 



where we have defined the dimensionless parameter 7 j/ETh, where Etyi = D/L^ is the Thouless 
energy on the scale L. Note that the total length of the wires which form the graph is Ltotai = '2'L{Nd + 1). 
Thus, the probabilities obtained via inversion of the matrix (10), cf. Eq.(9), are normalized on itotai- For 
further calculations, it is more convenient to change this normaUzation from Ltotai to the actual length of 
the system, 4L: Firstly, we recall that all Nd wires in the dot connecting the same two vertices have the 
same length, i.e. these wires are identical. Consider a point X inside the dot which belongs a given wire 
(out of Nd) and is infinitesimally close to one of the vertices a = 4 or 5. The probability to reach X from 
any other point is equal to the probability to reach a itself. Let us now introduce a probability V to reach 
X belonging to any of the Nd wires: 



(13) 



here N'^^^ ~ Nd if /3 is a vertex lying in the dot and N^^^^ = 1 otherwise. V is normalized on 4L and it 
reflects an enhancement of the probability for an electron to stay in the dot by the factor Nd- 

Furthermore, we define the piecewise continuous function V^{x,y) of continuous variables x,y G 
[— 2L, 2L] by selecting two appropriate vertices and replacing the wire-length parameters, Za, by x/L or 
y/L. E.g. the probability to reach any point y e [~L, L] in the dot from a point x E [~2L, —L] in the left 
contact, is given by V-y{x, y) — Ndj){M.'*)il after replacing Zi by x/L and Z4 by y/L. 

An analytic expression for V^{x, y) can be evaluated efficiently, but it is lengthy and will be published 
elsewhere. Besides, the inverse Laplace transform of 'P-^{x, y), cf. Eq.(4), can be calculated by exploiting 
the fact that all poles of V^{x, y) are simple and coincide with the zeros of the determinant of TW". Direct 
calculation yields^ 

(14) 



7), of the graph [11], implying 



detTW cx 5(7) = sinh(2v/^) {{Nd-l) + (A^d + l)cosh(2 V7) 



" We note in passing that S('y) is proportional to the so-called spectral determinant, det(— D A - 
that it does not depend on any of the auxiliary coordinates Zi . 
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Solving equation S{'^) = yields the following poles for the graph under consideration: 

7fe = -(^Y) 7fc = - (^kir + arccosy j;^^^ ^keZ. (15) 

Note that there is no pole at 7 = since the system is open. Defining the dimensionless function 

^(^^^^^^^l^,^,,^.^^ (16) 

where S'{x) — dxS{x), we can evaluate the time-dependent probability using the residue theorem by 
closing the integral contour in Eq.(4) on the left half-plane; 

V{x, y,t) = jY. ^(^' exp(7feSTht) . (17) 

k 

V{x, y, t) is plotted in Fig. 2 for fixed t = rTh/4, Nd = 2 and x either in the left contact or in the dot. We 
emphasize that for Nd > 1, 'P(-t, y, t) is discontinuous at ?; = ±L, describing confinement in the dot. In 
particular, V{x, y, t) = NdV{y, x, t) for a: in a contact and y in the dot. Normalization is reflected by the 
fact that V{x, y, t) satisfies a semi-group relation 

/2L 
dyVix, y, h)V{y, z, h) = V{x, z, h + h) . (18) 
-2L 

In the next sections we will evaluate the correction to the conductance and the dephasing time using the 
probabihty V. 



5 Quantum corrections to the conductance for the quantum dot model 

The classical conductance of the system described by Eq.(lO) is obtained via Kirchhoff 's circuit laws, since 
the contacts of length L and the central region of length 2L (with Nd wires in parallel) are connected in 
series. Denoting the contact conductance (i.e. the conductance of the left or right wire) as gc, we obtain 

5o = y(l + l/iVd)-i. (19) 

Note that the value of g^, cannot be chosen arbitrarily: Assuming that the substrate, from which the wire 
(length L and width W) is constructed, is 2D or 3D with mean free path £ and Fermi wavelength Xp, the 
conductance is given by g^^ = £W/ XpL or g'^^ ~ 21W^ /iTiX^pL. Our theory requires > 4/3 in order 
to obtain g > Ag, and quasi- ID diffusion requires Xp £,W L. For a quantum-dot of the size of 
several fim, etched on a GaAs/AlGaAs heterostructure (Xp O.OS/nm), we can estimate a typical value 
of 5c ^ 5. 
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To evaluate the quantum corrections Ag, Eq.(l), we need the return probability defined via Eq.(13) at 
coinciding a and (3. In this section, we consider the case T = (i.e., r^p oo) and study Ag as a 
function of the dissipation parameter 7. We calculate matrix elements [(A^''')^^)]ii, [(A^''')^^)]44 which 
yield the return probability for the dot: 

1 



ce[-L,L] 2VTD(iV<j + 1)^(7) 



(20) 



71 --2 



L 



and for the left wire: 

{^X ^ X^ 



xe\-2L,-L] 



X (iVd + l)sinh(v7(^ 
sinh(V7(f + 2)) 
X {Nd-l)(iNd + l)smh 



l)sinh 
(7Vd-l)sinh(V^|) 



(21) 



L 



-(7Vd-l)sinh| 

- (A^d + l)sinh(v/f (z+^)) ) " (^d + l)'sinh(v/f(|-2)) 

respectively. 7'-y(a;, a;) for the right wire, x G [L,2L], can be obtained from the symmetry property 
Vj{x, x) = V^i—x, —x). In the limit 7^0, Eqs.(20) and (21) reduce to 



'Po{x,x) 



L((iV, + l)2-(f)2) 



, 'Po{x,x] 



i(2 + f)(2-iV,f) 



(22) 



xel-2L.-L] 2D{Nd + l) 

00, since the central region is 



c6[-L,L] 2D{Nd + 1) 

Note that the return probability diverges for x G [—L, L] in the limit Nd 
effectively closed in this limit. 

Similarly to Eq.(19), the total quantum corrections have to be properly weighted by using the circuit 
laws. The total correction can be written as a sum over all wires i of the network [15]: 

1 v-^ dC 



A.9 



-4L> 



£2 dL, 



dxT^y{x, x) . 



(23) 



wire No. i 



where C is the effective total length of the system obtained similar to the total resistance. In the case under 
consideration, we have 

1 



C = L„ 



L 



2L{1 + 1/Nd) , 



(24) 



1/Li + • • • + 

where Lq = L corresponds to the left wire, LjVd+i — Lto the right wire and Li . . . L^^ ~ 2Lto the 
wires of the dot. We obtain the following expression for the total quantum correction: 



A.9 = -E- 



1 



Th- 



1 + l/Nd^ 



2L 



1 

dxV-y{x,x) + —2 I dxVj{x,x) + 



2L 

dxV^{x, x) 



(25) 



In Fig. 3, we show the total correction to the conductance according to Eq.(25) as a function of the dissi- 
pation parameter 7 for different values of Nd- Note that for 7 1 all curves are oc 1/^/^, similar to an 
infinite wire with different prefactors corresponding to different 

effective wire width. We are mainly interested in the regime 7^1, where the main result originates 
from the left and right wire and all curves approach the ergodic limit lim-y^o(A.g) = 2/3, cf. Eq.(2). 
This limit can be checked in this model by substituting Eqs.(22) into Eq.(25). Thus, in the absence of 
dissipation, our model has qualitatively the same behavior as RMT theory. Since the time to reach one 
contact from the other increases linearly with iV^, there is an intermediate regime at 1/Nd < 7 < 1 for 
large Nd, where the system is described effectively as two wires connected in series via the dot, which just 
plays the role of an additional lead, such that Ag = 1/3. 
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6 Evaluation of the dephasing time for the quantum dot model 

The dephasing time, r^, can be calculated from the phase difference acquired by an electron in a time- 
dependent (fluctuating) potential V{x, t) during a time-reversed traversal of its trajectory x{t) [9]: 

<S>[x{t)] = /" dr [V{x{t),t) - V{x{t), t ~ t)] . (26) 

"'0 

When averaged over the Gaussian fluctuations of the potential (e'*)\/ = e^^^* ^ , Eq.(26) leads to an 
exponential cutoff of the return probability^ 

V{x, X, t) r{x, X, t) ■ (e'^*[-(-)l){,(,)} « V{x, X, t) ■ e-^(--*) (27) 

where (. . means the average is over closed trajectories x{t) of duration t, staring and ending at 

X, and we defined the decay function F [16, 17]: 

F{x,t)= f dhJ{VV){x{t^),x{t2)M-t2)-{VV){x{h),x{t2),t-h-h)) .(28) 



^0 \ / {x(r)} 

In the case of the graph model for the quantum dot, the usual operational definition of reads 

F{x,T^{x))^l, (29) 

such that the correction to the conductance is given by Eq.(25) with a position dependent 7(2;) = 1/t^{x). 
The correlation function {VV) entering Eq.(28) is well known for the case of electron interactions in 
macroscopically homogeneous disordered systems [9]. Recently, we have generalized this theory for in- 
homogeneous, multiply-connected systems [10]. It has been shown that lyV) generically is given by 

47rT 

{VV){x,y,t)^—Po{x,y)5T{t), (30) 
where Po{x, x) ~ lim^^o P-i{x, x) and 

Srit) = nTwinTt) with w{x) = ^gQth(xJ - 1 ^^^^ 

sinh X 



•3 

Note that in the second equaUty of Eq.(27), we exponentiate the average over closed path, see Ref.[16] for details. 
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is a broadened (5-function which allows us to take into account the Pauli principle [16]. 
Inserting Eq.(30) into Eq.(28), we find 

AttT /■* 

J'ix, t) = / dh-2 Q{x, t^.tM - tm, t - Im) [Sriti - t2) - Srih + t2 - t)] , (32) 

9c Jo 

where tm = min[ti,t2] and = max[<i,i2]- The function Q is given by the dimensionless quantity 
DPq/L, averaged over closed random walks: 



J-2L V(Xo,Xo,ti+t2+t3) L 

All probabilities in Eq.(33) can be evaluated analytically from Eq.(17), Eq.(13) and Eq.(lO), by deriving 
the corresponding entries in the inverted vertex matrix [A^'']^^. The integrand is lengthy and we have 
chosen the following strategy for calculating the integrals; 1) We use Eq.(17) to rewrite Eq.(33) as: 

QixoM.h) = E n.f^^'^'^l^^l^^'^'I. , 6^-^^+^-^-+^'^- , with (34) 

n , A: , / 

-^7^(a;o, xi, 7l)7^(xl, X2, 72)7^(x2, xo, 73)-Po(a;i, 2:2) ■ (35) 

The integrals in Eq.(35) over space are evaluated symbolically with the help of a computer algebra program. 
2) Since the time dependence of Q in Eq.(34) is simply exponential, one of the time-integrals in Eq.(28) is 
calculated analytically. As a result, F{xo, t) simplifies to a single time integral and multiple sums: 

^ Qixo,jn,^k,li) / dr£(T,7„,7fc,7i) 

F(a;o, t) = — -^r— ■ (36) 

9c 2j7e(xo,xo,7m)e^'"^^''* 



V{xo, Xi,ti)V{xi,X2,t2)V{x2,Xo,t3) DPo{xi,X2) 



Here, the remaining time dependence of the kernel is incorporated in the function £: 



'S(r,7i,72,73) = u;(7rT)e'=i 



citT 



sinh{c2{tT-T))e-=3- _ sinh(ia(tT-r))cosh(e2T)e''3('r-^)/2 

2C2 C3 



, (37) 



with ci = (7i+73)fEf, C2 = (7i-73)|i^, C3=72^-Ci. (38) 

3) The sums and the integral over t are calculated numerically. 

This strategy allows us to calculate r^p and to describe the T-dependence of Ag in the quantum dot 
model, including the full crossover between different regimes of dephasing. 



7 Examples of application 

In this section, we use the graph model of the quantum dot to calculate r^p {xq , T) and Ag{T) in the case 
gc — 5 for the parameter Nd ranging from Nd — 1 (no confinement in the central region) to Nd — 100 
(almost closed quantum dot connected to ideal leads via two contacts). Our model is vaUd for this choice 
of gc, see the discussion in Section 5, and the total conductance of the system 1.25 < go < 2.5 is close to 
experimental setups [18, 19]. The results are shown in Fig. 4. 

The dephasing time is shown in Fig. 4(a) for several values of the origin of the Cooperon, xq, which can 
belong either to the central region (solid blue lines) or the contact (dashed red lines). To check the validity 
of the results, we compare r^p at high and small temperatures with earlier results for an almost isolated 
quasi- ID ring of total length 4L and total conductance gi [20]. 
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Fig. 4 (color online) (a) The dephasing time in units of the Thouless time, l/i^Th, plotted for several values of 
Nd and xo. Solid blue lines, correspond to xo — — 0.05L close to the center of the dot, while dashed red lines 
correspond to xo = —1.551/ close to the center of the left contact. The thickness (and brightness) of the curve 
determines the number of channels in the dot, Nd ~ 1, 10, 100 from thin to thick (and bright to dark). The black 
dotted lines correspond to the asymptotic results, Eqs.(39) and (40), derived from an isolated ring geometry, see main 
text for details, (b) The difference Ag + 2/3 between the correction to the conductance, Ag, and its universal zero- 
temperature value, Ag{T = 0) = —2/3, plotted as function of temperature. OD behavior of the dephasing time, 
characterized by Ag oc T^, appears at very low temperatures, requiring a precision much larger than 1% on the 
conductance measurement. Inset: (c) Total correction to the conductance —Ag (without subtracting Ag{T — 0), 
plotted as function of temperature. 



If tt Tip TTh = l/£'Th, dephasing is not sensitive to the boundary conditions and it is described 
by the theory of infinite systems [9]. In the ring, the high-T regime appears at T ^ giE^i^. The formula 
for Tip in this regime, including sub-leading terms, reads [20]: 



TTh V T^^^^T J 1 Stt'^^ ' ^' I 2gi / \ T J 90F VTri / V 



IC(1/2)I hr- +77% P^ • (39) 



We have reproduced this high-T behavior in the quantum dot model, see Fig. 4(a): Numerically obtained 
curves coincide with Eq.(39), after substituting NdQc for gi, when T ^ {gcNd)ETh- We note that 
dephasing in the high-T regime is substantially inhomogeneous in space, since the relevant trajectories are 
restricted to a small region around xq. In particular, for sufficiently high T, all curves for dephasing in the 
contact {Nd = 1, 10, 100) coincide with the curve for Nd = 1 in the central region, since the number of 
channels in the central region is irrelevant for dephasing in the contact. On the other hand, dephasing in 
the central region itself becomes weaker with increasing Nd, since Nd increases the effective conductance 
in this region. 
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In the low-T regime/ rxh ^ tt ^ Ti^, typical electron trajectories explore the whole system many 
times before dephasing becomes effective [1]. The geometry of the system is not important in this case 
and, therefore, the low-T regime is usually referred to as the regime ofOD dephasing. In the ring, it occurs 
at T < E'xh with given by [20] 

^ . HMl f f 1 + Jl + 1^ (40) 

TTh 327r2 \ T ) \ 45giETh 105 \EtiJ J 

The quantum dot model shows similar behavior at T ^ ^'xh- after substituting gc^d for the ring conduc- 
tance. We emphasize that OD dephasing in our model is governed by atypical trajectories, which explore 
the dot and the contacts many times during the time scale t ^ Tdw Therefore, the dephasing time is nearly 
coordinate independent: Dephasing in the central region and in the contacts is essentially the same. 

The correction to the conductance is shown in the inset. Fig. 4(c), for Nd = 1, 10, 100. We calculated 
Ag from the integral in Eq.(25) with a position dependent 7(2;) — 1/t^{x,T). As expected from the 
discussion in Section 5, the curves saturate to the universal value Ag ~ —2/3, when 7 = 7/-ETh ^ ^/Nd- 
Since 1/7 = t^/ttii ~ {gcNd){ETh/T)'^ in this regime and gc is small and fixed, saturation occurs when 
T < i?Th- The intermediate regime for 1 /Nd ^7^1, where Ag = — 1/3, cf. Fig. 3, is strongly distorted 
since it lies in the crossover region between high-T and low-T regime. We note that at T < lOi^xh, curves 
for different Nd look very similar Moreover, dephasing is very weak at T ^ i?Th where Ag is governed 
by a dwell time, Tdw, of the entire system and is practically T-independent. After subtracting the curve 
from its universal value, see Fig. 4(b), OD dephasing reveals itself as Ag cx for very low temperatures 
T < 0.2£'Th- At 0.2i?Th < T < i?Th one can observe only a transient, since (i) dephasing is not 
yet sufficiently weak to justify Ag oc and (ii) the OD regime of dephasing is not fully reached, cf. 
Fig.4(a). Moreover, if the leads are not perfectly absorbing, the transient can be extended even to lower 
temperatures due to additional dephasing in the leads. All this clearly shows that OD dephasing cannot be 
discovered directly in transport measurements through the quantum dot. Even at T < 0.2i?Th, a fitting of 
the experimental data would require g to be measured with a precision of much better than 1%. Alternative 
possibilities for the experimental observation of OD dephasing are discussed in the Conclusions. 

8 Conclusions 

We have suggested a graph model, which allows one to describe transport through mesoscopic quantum 
dots. The graph includes three quasi-lD regions: identical left- and right- wires and a central region. The 
identical wires are connected to ideally absorbing leads and mimic the contacts of a real system. The 
number of conducting channels in the central region can be of the order of- or substantially larger than the 
number of channels in the contacts. The latter case corresponds to a strong confinement of electrons in the 
central region. Thus the graph model is able to describe a crossover from opened to closed quantum dots. 

The model which we suggest can be viewed as complementary to the seminal RMT model. Firstly, the 
exact solution to the diffusion equation can be found for the graph model. Secondly, we have shown that 
our model correctly reproduces the universal regime of transport in full analogy with the RMT solution. 
Even more importantly, the graph model allows us to take into account interaction induced dephasing in a 
broad temperature range, i.e., we can describe the full crossover from ID to OD regimes. 

Using the solution to the diffusion equation on the graph, we have described in detail how to calculate 
the dephasing time and the weak localization correction to the conductance. Though the intermediate 
equations are rather lengthy, we have suggested an efficient combination of analytical steps (involving 
computer algebra) and numerical integration, which helped us to overcome technical difficulties. 

The general approach has been illustrated for the system with gc = 5. We have demonstrated that OD 
dephasing T^), which is governed by the Pauli principle and is very generic, occurs in the system at 



The intermediate regime, tt ^ T^h Tc/j, chai'acterized by oc T^^ is strongly distorted in the quantum dot, since: 
(a) The conductance gc is relatively small, reducing the range of validity of this regime, and (b) it occurs when typical electron 
trajectories are of the order of the system size making sensitive to the inhomogeneities of the graph. 
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T <C -Exh at arbitrary ratio of the channel numbers in the dots and leads. In this regime, dephasing is 
governed by atypical trajectories which explore the dot and the contacts many times during the time scale 
t ^ Tdw where the conductance is governed mainly by the dwell time and is almost T-independent. Our 
results confirm that weak OD dephasing is substantially distorted by the influence of the contacts and the 
leads. Therefore, its direct experimental observation in transport through the quantum dot would require 
not only very low temperatures but also unrealistically precise measurements. We conclude that alterna- 
tive experimental approaches are needed, where either the effects from the environment are reduced or 
the system is closed. One possibility to improve the effective precision of the measurements is related 
to extracting from the T-dependence of the Aronov-Altshuler-Spivak oscillations of the magnetocon- 
ductivity in almost closed mesoscopic rings. This option was discussed in recent papers [6, 20] where all 
effects of the environment were taken into account via a constant dwelling time. We plan to study in more 
detail the sensitivity of the AAS oscillations on the distortions from the environment using a ring model 
similar to the model of the dot presented here [21]. The other option is to extract r^p from experimental 
measurements of the electric or magnetic susceptibility of closed mesoscopic systems, e.g. by measuring 
the properties of resonators in which mesoscopic samples are deposited [22]. In closed systems, there is 
no universal limit of the quantum corrections at r^p ^ rd„, typical for transport through opened systems. 
Therefore, the saturation in the closed system can occur at much lower T, making them more suitable 
for an experimental observation of OD dephasing. A theoretical description of such experiments will be 
published elsewhere. 
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